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Excess emission over expected diffuse astrophysical backgrounds in the direction of the Galactic center region 
has been claimed at various wavelengths, from radio to gamma rays. Among particle models advocated to 
explain such observations, several invoke interactions between dark matter particles and ordinary matter, such 
as cosmic rays, interstellar gas or free electrons. Depending on the specific interstellar matter particles’ species 
and energy, such models predict distinct morphological features. In this study we make detailed predictions for 
the morphology of models where the relevant electromagnetic emission is proportional to the product of the dark 
matter density profile and the density of interstellar matter or cosmic rays. We compare the predicted latitudinal 
and longitudinal distributions with observations, and provide the associated set of relevant spatial templates. 


I. INTRODUCTION 

With the advent of large-scale sky surveys at frequencies 
spanning most of the electromagnetic spectrum, from radio 
to gamma rays, om' theoretical understanding of astrophysi¬ 
cal diffuse electromagnetic emission processes is confronting 
the test of observation in an unprecedented detailed way. Per¬ 
haps not surprisingly, at several frequencies diffuse emission 
models have, at times, fallen short of providing a satisfac¬ 
tory match to observations. Interestingly, in many cases such 
shortcomings are centered towards the inner regions of the 
Galaxy, a rich, and relatively poorly understood region. Over 
the years, some such excesses have found plausible explana¬ 
tions in the realm of “traditional” astrophysical processes, or 
in previously under-estimated emission from populations of 
astrophysical objects. In some cases, however, it has also been 
argued that the detected excess might have a “non-traditional” 
origin, possibly connected with new physics. 

Perhaps the longest-standing and most widely known such 
excess is the 511 keV line detected from a broad angular 
region by INTEGRAL/SPI HI- At larger energy, COMP- 
TEL reported an excess across the energy range between 1- 
20 MeV 0. Diffuse X-ray emission from the Galactic bulge 
region, with an approximately thermal spectrum with a very 
large associated plasma temperature (around 10 keV), has also 
been reported from Chandra data after point-source subtrac¬ 
tion El. At radio frequencies WMAP revealed excess mi¬ 
crowave emission at frequencies between 23 and 61 GHz, an 
excess known as WMAP haze 0. This radio “haze” has also 
been confirmed with Planck observations 0. Einally, sev¬ 
eral groups have identified an extended excess of gamma rays 
from the Galactic center region and from the inner Galaxy, in 
the few GeV range ||6H^. 

While astrophysical counterparts have been identified that 
might explain in part or entirely the excesses listed above, 
several studies have focused on the possible connection of 
the observed excess emissions with new physics, and specifi¬ 
cally with particle dark matter. Dark matter pair-annihilation 
or decay produces electromagnetic emission both as a result 
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of prompt photon emission from, e.g., neutral pion decay or 
internal bremsstrahlung or loop-mediated direct annihilation 
into photons, as well as from secondary mechanisms: the 
latter mechanism depends on how electrons and positrons, 
produced in dark matter annihilation or decay, loose energy 
via synchrotron, inverse Compton, Coulomb scattering and 
bremsstrahlung ISEl. 

As far as the primary emission is concerned, the predicted 
morphology follows the integral along the line of sight of the 
dark matter number density (squared) for decay (annihilation, 
respectively). Eor the secondary emission, instead, the mor¬ 
phology is complicated by the magnetic field structure and 
gas and electron densities ini. 

Alternately, some of the excesses listed above have been as¬ 
sociated with slightly less trivial new physics models, where 
the dark matter electromagnetic emission effectively depends 
on the environmental cosmic-ray population. Such models in¬ 
clude, for example, antiquark nuggets nmH, which would 
emit via several different mechanisms: free electrons would 
annihilate with positrons in the nugget’s electrosphere pro¬ 
ducing a 511 keV line; more energetic cosmic ray electrons 
would penetrate deeper and potentially produce photons in 
the COMPTEL energy range; cosmic-ray protons penetrat¬ 
ing into the quark matter would produce hadronic jets poten¬ 
tially responsible for Bremshstralung emission in the X-ray 
frequencies relevant for the Chandra excess; for proton cosmic 
rays penetrating deeper in the nugget, the complete absorption 
would eventually yield thermal photons with energies in the 
WMAP haze range (for a detailed review of all these mecha¬ 
nisms, see Ref. mi 

In eXcited dark matter (XDM) HU scenarios an excited 
state exists at energies of a few MeV above the ground state. 
Such state might be populated by collisions of the dark mat¬ 
ter particle with e.g. Galactic cosmic rays, yielding subse¬ 
quent electron-positron pairs that could explain the 511 keV 
line signal ua. 

Other scenarios where electromagnetic emission originates 
from elastic scattering of dark matter particles off of cosmic 
rays (or vice versa) were considered in Ref. USHnl, albeit 
no specific connection with any of the diffuse excesses was 
attempted. The relevant cosmic-ray populations are, in this 
case, high-energy cosmic-ray protons and electrons 113 . 

The generic feature of the class of models we focus on here 
is that the electromagnetic emission is proportional to the in- 
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tegral along the line of sight of the dark matter density times 
the density of a charged cosmic ray species. In this paper, 
we consider four classes of such cosmic ray species: (i) low- 
energy protons; (ii) low-energy free electrons; (iii) intermedi¬ 
ate energy (1 GeV) cosmic-ray electrons and protons; and (iv) 
high-energy (1 TeV) cosmic-ray electrons and protons. We 
then compare the predicted average longitudinal and latitudi¬ 
nal intensity profiles for the four cases (in some instances we 
even adopt more than one model for a given case) both with 
the morphological prediction for dark matter annihilation and 
with the observed excesse’ emission intensity. 


The remainder of this study is structured as follows: in 
sec. we describe the dark matter and cosmic-ray densities 
we employ in our analysis; sec. [^details on the calculation 
of the resulting emission profiles, while in sec. IV we com¬ 


pare our findings with the observed excesses’ profiles. The 
final sectionj^presents our discussion and conclusions. 
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In Figure [T] we show, in the right-most top panel, the inte¬ 
grated line-of-sight dark matter density squared (i.e. the mor¬ 
phology corresponding to an annihilating dark matter candi¬ 
date) for the case of the NFW profile. This morphology is of 
course azimuthally symmetric, and sharply centrally peaked. 
In the following two sections we will compare this benchmark 
scenario to more exotic morphologies arising from dark mat¬ 
ter interactions with cosmic-ray protons and electrons, inter¬ 
stellar gas, and free electrons. 


II. DENSITY DISTRIBUTIONS 


In this section we describe each of the matter distributions 
used throughout this study, including the density distribution 
of dark matter (sec. S’ as motivated by the results of cur¬ 
rent generation N-body simulations; the density of Galactic 
cosmic-rays for energies between a few hundred MeV and 
several TeV (sec. II Bi, as derived through numerical simu¬ 
lations of cosmic-ray propagation; the distributions of inter¬ 
stellar gas (sec. m in two recent and distinct modeling ap¬ 
proaches; and, finally, the most up-to-date model of the Galac¬ 
tic distribution of free electrons (sec. m- For dark matter, 
cosmic rays, and gas we also discuss the dominant sources of 
uncertainty and attempt to bracket the range of state-of-the art 
models. 


A. Dark Matter 


We employ as our benchmark dark matter profile a Navarro- 
Frenk-White (NFW) MU profile with inner slope a = 1 and 
scaling radius = 20 kpc. In order to assess the uncertain¬ 
ties due to choice of dark matter profiles, we also utilize a 
steeper (more ‘cusped’) generalized NFW with an inner slope 
of a = 1.2; this is in part motivated by recent measurements 
of the GeV excess M; a steeper inner slope is physically 
motivated in the context of halo evolution including adiabatic 
contraction. On the opposite extreme, we consider a cored 
Einasto ||T91 profile with aE — 0.16. We are not concerned 
with relative normalizations here since we are exclusively in¬ 
terested in the morphological predictions of the models under 
consideration, not the overall emission intensity. The func¬ 
tional forms of the two dark matter density profiles we con¬ 
sider here are as follows: 


B. Interstellar Gas 

Some models of dark matter - see e.g. Ref. Cl for a re¬ 
view - predict interactions with Galactic gas and/or free elec¬ 
trons. Both of these distributions are strongly peaked toward 
the Galactic center (GC), and a detailed understanding of the 
inner few kpc of the Galaxy are needed to formulate solid pre¬ 
dictions for the resulting morphology. 

The gas density in the Milky Way is typically described by 
summing contributions from three dynamically distinct com¬ 
ponents of hydrogen gas in molecular, atomic, and ionized 
phases. The former two overwhelmingly dominate the gas 
density near the GC, and are therefore the most important 
here. 

The three-dimensional distribution of gas in the Galaxy 
can be determined with excellent accuracy by combining sur¬ 
veys of atomic transition lines with a Galactic rotation curve. 
For atomic hydrogen, the hyperfine transition at 21 cm pro¬ 
vides a direct observable. In the optically thin limit the col¬ 
umn density is related to the observed brightness tempera¬ 
ture by a single parameter: the hydrogen spin temperature 
Ts- In the case of molecular hydrogen, the lack of a perma¬ 
nent dipole moment requires, instead, use of a tracer gas, the 
CO( J = 1 —> 0) transition, which is related to the molecu¬ 
lar hydrogen density through a conversion factor Xco- This 
factor is, in principle, spatially varying. 

The deconvolution technique described above relies on a 
relative velocity between the gas and the solar system. In the 
direction of the GC, the gas is co-rotating, implying no kine¬ 
matic resolution, leading to a distance degeneracy along lines 
of sight near Galactic longitudes I « 0. This problem is com¬ 
pounded by the so-called ‘near-far ambiguity’ which corre¬ 
sponds two distances to the same radial velocity in the inner 
Galaxy. In order to alleviate such problems, one can incorpo¬ 
rate a model into the deconvolution procedure. For CO (H 2 by 
proxy), we use the gas model from Pohl, Englmaier and Bis- 
santz (PEB), Ref. Il20l . which combines the survey of Dame 
et al ETIl with a gas flow model derived from hydrodynamic 
simulations as well as interpolation of the spiral arms across 
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the line-of-sight toward the GC. Not only does this provide a 
significant improvement in the l.o.s. gas distribution toward 
the Galactic center, but it importantly reconstructs a promi¬ 
nent Galactic bar and the intervening spiral arms. For atomic 
hydrogen, Nakanishi and Sofue (NS), Ref. Il22l . assume, in¬ 
stead an analytic model of hydrodynamic equilibrium for the 
scale height of HI as a function of Galactic radius and uses 
the Leiden-Argentine-Bonn survey ll^ 21 cm survey. This 
combination provides a high resolution model of H 2 and HI 
which we denote PEBh-NS. 

One can alternatively attempt to build a three-dimensional 
model based on observations of individual structures that are 
prominent in the region. This is the approach of Refs. oca 
which predicate two primary disk components in the inner 
3 kpc of the Galaxy: a dense “Central Molecular Zone” ap¬ 
proximately centered on the GC and aligned with the Galac¬ 
tic Plane, as well as a holed “Galactic-Bulge-disk” which is 
rotated at 13.5° counter-clockwise to the Galactic plane as 
well as having it’s major axis inclined « 45° away from the 
line of sight, which is notably much larger than the PEBh-NS 
case, leading to a larger projected extent. Beyond the in¬ 
ner 3 kpc, we use an azimuthally averaged gas profile from 
Ref. Il26ll to describe the molecular phase, and the HI pro¬ 
file from Ref fTf\ . corrected to the updated normalizations of 
Ref ll28l . We collectively refer to this model as Eerriere 2007 
(E07). 

Eor both the PEBh-NS and E07 models, we incorporate the 
contributions of dark gas - i.e. molecular hydrogen not traced 
by CO emission, and atomic hydrogen missed due to the as¬ 
sumption of a uniform spin temperature - by re-normalizing 
the total ‘analytic’ HI column density along each line of sight 
to that of the GALPROP map Essentially, this involves 
a fitting, in addition to HI and H2, a dust template 1291 to 
gamma-ray data. The detailed construction of this template is 
described in detail in Refs. Il30l[3l]. 

Einally, for both models, we also include a contribution 
from warm, hot, and very hot phases of ionized hydrogen 
based on the NE2001 l3^ 1^ model for free electrons, de¬ 
scribed below. The specific implementation we use is de¬ 
scribed in Refs. Esiia. 

We note that the precise gas distribution in the outer galaxy 
is of only marginal importance due to the central peak in the 
dark matter (DM) density profile. Consequently, we assume 
both Ts = 150if and Xqo = 2 x cm“^(K km/s)“^ 
and neglect spatial dependence of either one of these quanti¬ 
ties. Although spatial variations are certainly present in both 
components l34l[T5l . we expect the effect of this on the mor¬ 
phology of the Gas x DM profile to be small in the GC region 
of interest. 


C. Cosmic-Rays 

In order to obtain the steady state distribution of cosmic- 
rays in the Galaxy, we use the numerical code GALPROP 


* rbandsjii 12_v2_qdeg_zmax 1 _Ts 150_EBV_mag5_limit.fits.gz 


v54 .r2504|^ lf3^ fJTl , which encompasses all of the physics 
relevant for cosmic-ray transport through the galaxy in¬ 
cluding energy losses, primary source distributions, and re¬ 
acceleration. A detailed desertion of the physics can be 
found at the dedicated Web Sit^B 

Our default diffusion setup consists the model 
^1^4^20^150^5 in Ref fill. This model features a 
standard set of diffusion parameters fit to local observations 
of a variety of primary and secondary species including 
protons, helium, electrons, positrons, B/C, and ^°Be/®Be 
El- The simulations assume cylindrical geometry with 
free escape boundary conditions from a diffusion halo of 
radius 20 kpc and half-height 4 kpc. Because the dark matter 
profiles are strongly centrally peaked, most of the reasonable 
variations related to the diffusion setup and geometry are 
unimportant to the final cosmic ray (CR) morphology. This 
includes the halo height, provided z^ax is greater than a few 
kpc and the CR energy of interest is less than a few TeV, 
at which point CR propagation transitions to the rectilinear 
regime. Notably, since we are concerned with the morphology 
and not the spectrum of cosmic rays, parameters describing 
diffusive re-acceleration and injection spectra are irrelevant. 
The energy dependence of the diffusion constant will have 
only a very minimal effect on morphology at the very lowest 
and highest energies. Eor definiteness, we use here a diffusion 
coefficient L»(7^) = 5.3 x \Q'^^cin^s-^(n/A 

Of particular interest to this analysis is the distribution of 
primary cosmic-ray sources. In what follows, we show predic¬ 
tions for four tracers of supernova remnants, believed to gen¬ 
erate the majority of Galactic cosmic-rays. These are based 
on the observed surface densities of Galactic pulsars (Yusi- 
fov and Lorimer ||39l), OB stars m, and supernova rem¬ 
nants (SNR ||4T|). Each of these distributions suffer from sub¬ 
stantial uncertainties in the inner 1-2 kpc due to both statistics 
and systematics surrounding correction for selection effects. 
In fact, the latter three distributions are parametrized such that 
the surface density is zero at the Galactic center, while the 
Yusifov case is non-zero, leading to a potentially large differ¬ 
ence in CR densities in this region. Simulations of neutron star 
populations in the Milky Way’s gravitational potential also in¬ 
dicate a non-zero central density 1421 . with certain models 
predicting a strong peak. Einally, it is difficult to rule out 
the possibility of a cataclysmic event, or enhanced injection 
of cosmic-rays at the GC, as may be evidenced by the Eermi 
Bubbles ||43]| . and perhaps the Galactic Center 7-ray Excess 
itself (see e.g. Refs. 1B1I451). Eortunately, only for low ki¬ 
netic energies - and only for very high energies for the case 
of - does this have a significant impact on the CR distri¬ 
butions, which become substantially smoothed after diffusion. 
In what follows, we do not include these more speculative CR 
source distributions. 


^ Cun'ent versions of GALPROP are available at 
http://galprop.sourceforge.net/ 

^ http://galprop.stanford.edu 
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D. Free Electrons 

The density of free electrons Ue in the Milky Way has been 
mapped to reasonably good precision by fitting complex mul¬ 
ticomponent models to a thousands of measurements includ¬ 
ing primarily pulsar dispersion measures (the line-of-sight in¬ 
tegral of the rie), temporal and angular broadening of radio 
pulses sensitive to variations in rig, and emission measures 
(the line-of-sight integral of rig). Combining this huge body 
of information has lead to the development of the so-called 
NE2001 model 

NE2001 contains 5 component classes: (i) smooth compo¬ 
nents including a thin and thick disc as well as 4 logarithmic 
spiral arms, (ii) a Galactic center region consisting of a Gaus¬ 
sian component in radius and scale height, (iii) a 4 component 
local ISM region, (iv) 78 clumped HII regions, and (v) 16 void 
regions. The only region with high DM density is the Galactic 
center, and thus we only include contributions from the dense 
and thin (140 pc) disk, the lower density thick disk, and the 
small, but very dense Galactic center components. We neglect 
the spiral arms, which do not extend to the inner 3 kpc of the 
Galaxy as well as clumps and voids which provide only weak 
and small angular scale (6 ^ 1°) perturbations to the overall 
density profile in the GC direction. 

We also utilize a correction to the thick disk scale-height 
proposed in Ref. 1461 after recalibrating the NE2001 model 
while avoiding strong HII regions in the Galactic plane. This 
correction roughly doubles the height of the thick disk to 
« 2 kpc, which for our purposes leads to a substantial broad¬ 
ening of the expected emission to high latitudes and is there¬ 
fore less disk-like and more spherically symmetric. Ref. ia 
provides a recent comparison and recalibration of all models 
of the free electron density. Importantly, they are all based on 
approximately the same ingredients, with the primary differ¬ 
ences in the smooth components being the thick disk scale- 
height. The model used here is the thickest, and thus offers 
the most optimistic scenario for fitting a DM x free electron 
signal to an approximately spherical excess. As we will see 
below, the primary motivation for modeling such emission is 
to explain the sharp excess observed at 511 keV by INTE- 
GRAL/SPI and even this optimistic model is too disk-like to 
explain the signal. 

As a final note, the next generation of free electron mod¬ 
els are likely to incorporate all-sky surveys of the hydrogen-a 
emission line, perhaps leading to a significant modification of 
the model presented here. Under the assumption that the num¬ 
ber density of ionized hydrogen (i.e. free protons) is equal 
to the number density of free electrons, the Ha intensity is 
directly proportional to the emission measure ll48ll . Unfor¬ 
tunately, such corrections are non-trivial and prone to large 
errors due to the quadratic dependence of the integrand on rie. 


III. MORPHOLOGICAL PROLILES 

With the matter distributions as specified above, we can 
now perform the line-of-sight convolution of dark matter with 
our gas models (meaning the integral along the line of sigh 


of dark matter times relevant cosmic ray or gas density). In 
the top row of Eigurel^we show the projected emission pro¬ 
files for an NEW profile convolved with the E07 and PEBh-NS 
gas models, the NE2001 free electron model, and a standard 
NEW profile in order to compare against a standard annihi¬ 
lation morphology. In each case, use of either a contracted 
NEW or Einasto profile would lead to a slightly brighter and 
cuspier, or fainter and more cored profile in the innermost re¬ 
gions and each template has been normalized to its maximal 
value, with no other salient differences. 

The E07 gas model possesses a very bright central core due 
to the highly concentrated CMZ zone surrounding the galac¬ 
tic nucleus (not to be confused with the even denser circum- 
nuclearring ll49ll which occupies the innermost wlO pc and is 
below the scale probed here). Emission from the bulge disk 
can be seen in the diagonally oriented flares on either side. 
Beyond 5 degrees from the GC, the thin molecular and atomic 
disks dominate the emission and do not significantly extent to 
high latitudes making it nearly impossible to obtain a spheri¬ 
cal excess. Similarly, the empirically derived PEBh-NS model, 
shows a bright central disk which more smoothly falls off into 
the broader Galactic disks. In addition, the high latitude emis¬ 
sion is significantly enhanced with respect to E07. Still, the 
overall profile lacks azimuthal symmetry and we can conclude 
that a truly spherical excess is not well fit by DM x gas pro¬ 
file. 

The case of free electrons is more subtle. Here, we observe 
a thick disk extending beyond 2 kpc (b > 15°) which, com¬ 
bined with the DM halo produces a roughly spherical emis¬ 
sion profile. A very bright emission disk can be observed at 
the Galactic center, though the angular extent is less than The 
thin disk, however adds a distinct elongation along the plane 
making the averaged longitude profile significantly less steep. 

In the bottom row of Eigure[T]we show projected emission 
templates for benchmark cosmic-ray protons and electrons at 
1 GeV and 1 TeV, using a Yusifov profile for the primary 
source distribution. Eor low and high energy protons, as well 
as low energy electrons, the CR density is relatively uniform 
over the region of interest, leading approximately to the same 
profile as would be expected from dark matter decays. Eor 
electrons and positrons, inverse-Compton and synchrotron en¬ 
ergy loss timescales (Tics,sync oc E~^) limit the diffusion ra¬ 
dius to Rdis = \/D{E)Tics,sync « T.hkpc/EceV for our 
choice of parameters, implying that the CR distributions will 
depend significantly on energy above a few tens of GeV as 
the cosmic-rays lose energy before propagating farther than a 
few kpc from their production region (provided the primary 
source distribution is not uniform). Supernova remnants are 
highly concentrated in the plane of the Galaxy, resulting in a 
full-width-half-max of only a few hundred parsecs near 1 TeV. 
Therefore any dark matter model which predicts emission due 
to interactions with high energy electrons or positrons will re¬ 
sult in a significantly disk-like morphology. 

An important caveat should be kept in mind before exclud¬ 
ing models based solely on a disk-like component. Namely 
that most “excess” signals rely on fitting and subtracting off a 
complicated background model especially at low Galactic lat¬ 
itudes. In some cases, particularly those with gas correlated 
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FIG. 1. Top, from left to right: Projected emission profiles for DM (NFW profile) times the F07 gas model, PEB+NS gas model, free 
electrons, and, for reference a NFW^ profile (DM annihilation). Bottom: Projected profiles of DM times cosmic-ray protons and electrons for 
representative energies of 1 GeV (two left-most panels) and 1 TeV (right-most two panels) as calculated with the GALPROP package, using 
the Yusifov 1381 distribution of primary sources. 


backgrounds, this fitting procedure can potentially also sub¬ 
tract off a disk component which was actually part of the sig¬ 
nal. For example, the background models used to obtain the 
Galactic Center excess rely on fitting independent gas annuli 
to gamma-ray data without including any model for excess 
emission. Thus if the Galactic center excess actually contains 
a disk-like (gas correlated) component this could be ‘hiding’ 
in the artificially enhanced background normalization. One 
way to alleviate such issues is to fit all components simultane¬ 
ously rather than relying on residual emission alone. 

Finally, we provide skymaps in FITS file format as well 
as Python scripts at the supplemental materials web pag^ 
for each combination of halo profile/gas modePfree electron 
model, as well as logarithmically spaced templates in CR en¬ 
ergy. Details about the files and model assumptions are pro¬ 
vided on the same web page. We also provide example scripts 
demonstrating how to specify and integrate an arbitrary three- 
dimensional distribution against the cosmic-ray, free e~, and 
gas models which may be found useful in other contexts. 


'* http://planck.ucsc.edu/dmcr-morphology 


IV. COMPARISON TO OBSERVED GALACTIC CENTER 
EXCESSES 

We now proceed to comparing the morphology resulting 
from the models under consideration here to three residual 
signals observed in the Galactic center region. 

The first signal of interest is the 511 keV line. First detected 
in 1972 ll50ll . the 511 keV excess has since been the subject of 
lively scientific debate. The most reliable measurement m 
is provided by the SPI instrument on-board the ESA’s INTE¬ 
GRAL satellite which indicates a large excess population of 
low-energy positrons in the Galactic center region. Proposed 
explanations range from cosmic-ray interactions with the in¬ 
terstellar medium to exotic astrophysical objects, to dark mat- 
teij^ However, any plausible explanation must not only ex¬ 
plain the spectrum, but also the basic morphological proper¬ 
ties of the observed excess which are well well fit by two com- 


^ For an extensive review of this signal and its possible origins we refer the 
reader to Ref. ED 
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FIG. 2. Profiles in Galactic longitude (left column) and latitude (right column) of the 511 keV excess detected by INTEGRAL/SPI taken from 
Ref □. From the top to bottom row, we show projected profiles for the dark matter density convolved with various distributions (see text) of 
Galactic gas, free electrons, cosmic-ray electrons, and cosmic ray protons. The shaded regions show variations due to the choice of dark matter 
density profile for the benchamark NFW model (solid blue). Each profile has been convolved with a Gaussian with FWHM 3° to approximate 
the INTEGRAL/SPI point spread function, and has been normalized to the same value atl,b = 0. The profiles are averaged over |6|, |(| < 30°, 
respectively; notice that no error bars are given in Ref. (T|. 


ponents: (i) a spherical-bulge consisting of two radial Gaus- 
sians (projected FWHMs of 3° and 11°) and a r:! 30% fainter 
thick disk (FWHM 7°) 13^ . 

In Figure 1^ we compare longitude and latitude profiles m 
of the 511 keV excess in the Galactic center region against 
each of the four morphological profiles. To account for the 
finite angular resolution of the instrument, we have applied a 
Gaussian point spread function with a FWHM of 3°. The lon¬ 
gitude and latitude profiles have been averaged over |6| ,1^1 < 
30°, respectively, in order to match the same angular averages 
as in the reported excess of Ref. HI. Each distribution has 
been normalized independently in both latitude and longitude 
at 6 = 0 which is suitable for observations with a large point 
spread function. In each panel we also show the prediction for 
a squared NFW profile, for comparison to the standard case of 
annihilating dark matter. To represent variations on the dark 
matter halo, we show the canonical convolution with an NFW 
profile in dark blue, with shaded cyan (red) regions represent¬ 
ing the contracted NFW (Einasto) profiles. 


In the top row of Eig. |^we show the DM x Gas for our 
PEBh-NS and E07 models. As observed in the previous sec¬ 
tion, the latitude profile is extremely centrally peaked and is 
much to steep to match the observed excess. Likewise, the 
longitude profile is also very centrally peaked due to the very 
bright CMZ region surrounding the GC. This is particularly 
true in the E07 model, where off-center emission is orders of 
magnitude lower beyond 5°. Variations in the halo profile are 
relatively small and unable to reconcile a compatible signal. 

In the second row, we show the free-electron case. Here 
the halo model has a much larger impact due to the ex¬ 
tended distribution of free e~. Because 511 keV corre¬ 
sponds to the electron/positron mass, DMxEree is per¬ 
haps the primary case of interest for this signal, as DM mod¬ 
els producing an abundance of low energy positrons will form 
para-positronium bound states which subsequently decay into 
two 511 MeV 7 -rays. Another possibility arises in models 
of Quark ‘(anti)-Nugget’ Dark Matter [121, in which heavy 
‘nuggets’ of quark matter form as a result of a new high- 
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FIG. 3. Same as Figure]^ but for the 3-10 MeV COMPTEL excess of Ref. (3, after subtracting the average ‘base’ emission. Predicted profiles 
have been convoluted with a PSF of FWHM=3.8°, and have been averaged over |6| < 5° and 1° < |1| < 30°, respectively. 


density color-superconducting QCD phase. Such models can 
explain both baryogenisis and provide a connection to the 
« 5 : 1 relative abundance of dark matter and baryons. In 
these models quark matter can be surrounded by a layer of 
positrons which can annihilate through collisions with free 
electrons. As noted in the Introduction, these models also pre¬ 
dict a variety of other multi-wavelength signatures which are 
tested below and act as an interesting test case throughout for 
confirming or ruling out models based on the detailed mor¬ 
phology of DM-matter interactions. Our results indicate that 
for any angular scale larger than a few degrees, the predicted 
morphology for DM x Free electrons does not match the ob¬ 
served signal, with significant over-predicted intensity around 
10 degrees in longitude. 

In the third row, we show the emission profile for cosmic- 
ray electrons at 1, 32, and 10^ GeV and protons at 1 GeV 
and 1 TeV, along with the variations introduced by CR source 
distributions. In eXcited DM (XDM) lfT4l scenarios it has 
also been proposed that a DM particle is collisionally excited 
through some inelastic interaction with, e.g. cosmic-rays or 
gas, into a higher energy state. If the difference in energy lev¬ 
els is 5E = 2TOe, such a state could decay through emission 


of a non-relativistic e+e” pair, where the subsequently an¬ 
nihilates. Our results indicate that for either protons or elec¬ 
trons the predicted morphology is qualitatively incompatible 
in both latitude and longitude with the observed 511 keV sig¬ 
nal profile. The morphology predicted for pair-annihilating 
DM models provides a much better fit to the observed mor¬ 
phology. 

We now consider observations from the Compton Tele¬ 
scope (COMPTEL), which detects a significant and broad 
spectrum excess from 1-20 MeV which is well over the 
expected inverse-Compton and bremsstrahlung backgrounds. 
Limited statistics and systematics surrounding the deconvolu¬ 
tion techniques lead to a somewhat noisy morphological pro¬ 
file which is difficult to definitively call diffuse, as opposed to 
some combination of point sources (e.g. hard X-rays from su- 
pemovae). Models exist which postulate a sharp steepening in 
the cosmic-ray electron spectrum below 200 MeV exist which 
do not violate synchrotron constraints ll53l and can potentially 
also link the hard X-ray excess from INTEGRAL to the low 
energy 7 -ray regime (EGRET). However, such models lack a 
straightforward physical interpretation and are remain ad hoc, 
leaving open the door for new physics explanations. A variety 
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of astrophysical models are reviewed in Ref. 

Figure shows our results for the COMPTEL 3-10 MeV 
excess, where the average “base” emission has been sub¬ 
tracted off El. The predicted emission profiles are again nor¬ 
malized to the l,b = 0 value, and have been averaged over 
\b\ < 5° and 1° < |Z| < 30°, respectively and smoother 
out with a Gaussian PSF of FWHM=3.8° as warranted for 
the COMPTEL observations. The four rows and two columns 


adopt the same conventions as in Fig. IV 


In the context of the anti-quark nugget model, the COMP¬ 
TEL emission might be associated with non-relativistic cos¬ 
mic ray electrons penetrating deeper in the nugget’s electro¬ 
sphere. As discussed in Ref. El the relevant electron veloci¬ 
ties to explain the COMPTEL excess spectrum range between 
v/c = 0.001 and 0.01. Therefore, the relevant population 
corresponds to the free electrons case (second row), although 
it is interesting to also compare with the cosmic-ray electron 
case (third row). While the CR electron longitudinal prohle 
is marginally compatible with the flat COMPTEL profile, the 
free electron case is significantly under-abundant. The latitu¬ 
dinal profile, instead, shows a remarkably close match for the 
free electron case to the observed profile, but only in the inner¬ 
most few degrees, while falling short at larger latitudes. The 
cosmic-ray electron (or proton) latitudinal profiles are instead 
significantly different from the observed excess. Finally, con¬ 
volution with gas density profiles also provides a decent match 
to observations as a function of latitude in the innermost few 
degrees. 

In summary, models invoking a convolution of free elec¬ 
trons and dark matter to explain the COMPTEL excess, such 
as antiquark nuggets, provide an adequate morphological 
shape for the innermost regions at low latitude, but fall short 
at larger latitudes and provide a relatively poor match to the 
longitudinal profile. 

Finally, we come to the observed Fermi-LAT GeV ex¬ 
cess. Detected initially in 2009 1551 . a recent, detailed re- 
evaulation || 6 | has led to renewed interest and subsequent de¬ 
tection by several other groups |[8l|56l. Based on current state- 
of-the art diffuse 7 -ray background models, the excess is ap¬ 
proximately spherical with a radial profile consistent with a 
contracted NEW. The spectrum is also well fit by a « 50 GeV 
weakly interacting massive particle annihilating to 5-quarks 
with a cross-section near the thermal relic value, making the 
signal a prime candidate for the potential detection of dark 
matter. Thus far, proposed astrophysical explanations in terms 
of unresolved pulsars l57ll have failed to explain the hard low- 
energy spectum, luminosity, and spatial extent ll58l . which has 
now been detected beyond 15° ll56ll (see however, the recent 
Ref II 59 I which proposes a new population of very young pul¬ 
sars). A viable interpretation involving cosmic-ray outbursts 
from the Galactic center region MM remains an open ques¬ 
tion, though for hadronic models the significant gas correlated 
emission from neutral pion decay imposes substantial con¬ 
straints. 

In Figurejl^we exhibit our results for the latitudinal profile 
from the Galactic center taken from Ref. Q. Here, we smooth 
the predicted profiles with a PSF of FWHM=0.3°, and we av¬ 
erage, as in Ref. ||2l over |/| < 2°. Unlike the previous cases 



FIG. 4. Same as Figure]^ but for the ~GeV Fermi-LAT ‘excess’ 
from Ref. (7). Predicted profiles use a PSF of FWHM=0.3°, are 
averaged over |1| < 2°, and are normalized at 5 = -|-6°. 


we normalize the predicted emission profiles at 5 = 4 - 6 °. 


Our key findings are that models predicating a convolution 
of DM with interstellar gas or free electrons produce a remark¬ 
ably good fit to the latitudinal distribution of the observed pro¬ 
file (top two rows). Models, instead, involving DMxCR con¬ 
volutions do not provide a good match, over predicting sig¬ 
nals at large latitudes and falling short in the inner regions — 
in other words producing unacceptably shallow density pro¬ 
files. This rules out for example interpretation based on the 
models discussed in Ref. M- Although longitudinal distribu¬ 
tions are currently unavailable, the observed excess has been 
found to be highly spherical making such that any profile with 
a significant disk component is unlikely to be viable barring 
substantial changes to the diffuse background model. 
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V. DISCUSSION AND CONCLUSIONS 

We carried out a model-independent study of scenarios 
where diffuse electromagnetic emission originates from in¬ 
teractions of dark matter particles with the interstellar gas, 
free electrons or Galactic cosmic rays. We assumed that the 
relevant morphology depends on the integrated line-of-sight 
product of the dark matter particle density times the relevant 
gas/cosmic-ray density. The key motivation to consider such 
models is that a variety of large-scale diffuse excesses in the 
general direction of the Galactic center have been identified, 
at wavelengths ranging from radio to gamma rays, and that a 
variety of particle physics scenarios have been proposed that 
rely on dark matter interacting with Galactic gas or cosmic 
rays. 

We considered a variety of state-of-the-art gas density mod¬ 
els, and well-motivated cosmic ray models, utilizing several 
different assumed injection source distribution profiles. We 
produced the relevant latitudinal and longitudinal profiles, and 
we are making our results available publicly on the web. Fi¬ 
nally, we compared our predictions with the emission profiles 


of a few selected diffuse excesses. 

Generally, we find that the morphology of the 511 keV 
line is matched much more closely by dark matter pair- 
annihilation than by any of the dark matter times gas/cosmic 
ray models we investigated; the COMPTEL excess is repro¬ 
duced by models of dark matter times gas or free electrons in 
the innermost regions, at low Galactic latitude, but not at large 
latitudes or in longitude. Finally, models advocating dark 
matter/cosmic-ray convolution fail at reproducing the Fermi- 
LAT GeV excess, while dark matter times gas or free electrons 
provides a relatively good fit to the latitudinal distribution of 
the signal. 
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